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Abstract 

We report on an improved study of 13^ /)(*)o^(*)o decays, based on 274 x 10^ BB events 
collected with the Belle detector at KEKB. The following branching fractions have been measured: 
B{B^ D^K°) = (3.72 ± 0.65 ± 0.37) x 10^^ and B{B^ D^K*°) = (3.08 ± 0.56 ± 0.31) x 10"^ 
We also obtain evidence with 3.2c7 significance for S° D*^K^ with B{B^ D*^K^) = {?,.2t\l± 
0.4) x 10"^ No significant signal has been found for the S° ^ D*^K*^ and B° ^ L)(*)Oi^*o decay 
modes, and upper limits at 90% CL are presented. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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Within the Standard Model, CP violation arises due to a single phase in the Cabibbo- 
Kobayashi-Maskawa quark mixing matrix fH]. Measurements of the Unitary Triangle angles 
therefore test the consistency of the Standard Model. Precise measurements of the branching 
fractions for 5° ^ D^K*^, 5° D'^K*^ and B° D^^p^*^ decays, where D^^p denotes 
or decay to a CP eigenstate, will allow a measurement of the angle 03 0. The decay 
B'^ D^K^ can also be used to measure time-dependent CP asymmetry in B decays 3|. 
The 5° Di*)oxo g^^id B^ /}(*)o^*o decays have been previously observed by Belle jij. 

For all these measurements, knowledge of the ratio of amplitudes Rdor*" = a[bo^dok*o] ~ 
b[so^doR*o^) crucial importance; its value is predicted to be about 0.4 
Here we report improved measurements of B^ /}(*)o^o and search for B^ — >■ Z)(*)*'A'*° 
and S° D^*)^K*^ decays with the Belle detector at the KEKB asymmetric energy 
e^e~ collider 0. The results are based on a data sample collected at the center-of-mass 
(CM) energy of the T(4S') resonance, which contains 274 x 10^ produced BB pairs. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOP), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere f^. Two different inner detector 
configurations were used. For the first sample of 152 million BB pairs, a 2.0 cm radius 
beampipe and a 3-layer silicon vertex detector were used; for the latter 122 million BB pairs, 
a 1.5 cm radius beampipe, a 4-layer silicon detector and a small-cell inner drift chamber were 
usedQ. 

Charged tracks are selected with a set of requirements based on the average hit resid- 
ual and impact parameter relative to the interaction point (IP). We also require that the 
transverse momentum of the tracks be greater than 0.1 GeV/c in order to reduce the low 
momentum combinatorial background. 

For charged particle identification (PID), the combined information from specific ioniza- 
tion in the central drift chamber (dE/dx), time-of-fiight scintillation counters (TOP) and 
aerogel Cerenkov counters (ACC) is used. At large momenta (> 2.5 GeV/c) only the ACC 
and dE/dx are used. Charged kaons are selected with PID criteria that have an efficiency of 
88%, a pion misidentification probability of 8%, and negligible contamination from protons. 
All charged tracks having PID consistent with the pion hypothesis that are not identified as 
electrons are considered as pion candidates. 

Neutral kaons are reconstructed via the decay Kg —>■ vr+vr" with no PID requirements for 
these pions. The two-pion invariant mass is required to be within 6 MeV/c^ (~ 2.5a) of the 
nominal mass and the displacement of the tt'^tc" vertex from the IP in the transverse 
(r — (p) plane is required to be between 0.2 cm and 20 cm. The direction from the IP to the 
TT+TT^ vertex is required to agree within 0.2 radians in the r — (p plane with the combined 
momentum of the two pions. 

A pair of calorimeter showers not associated with charged tracks, with an invariant mass 
within 15 MeV/c^ (~ 3a) of the nominal 7r° mass is considered as a tt" candidate. An energy 
deposition of at least 30 MeV and a photon-like shape are required for each shower. 

K*^ candidates are reconstructed from K~it~^ pairs with an invariant mass within 
50 MeV/c^ (IP) of the nominal K*^ mass. The polarization of K*^ mesons in B decays 



4 



is also utilized to reject background through the use of the helicity angle 6k*, defined as the 
angle between the K*^ momentum in the B meson rest frame and the K~ momentum in 
the rest frame. We require |cos(6'x*)l > 0-3 for 5° ]ji*)ox*o decay channel. This 
rejects about 20% of the background and retains 97% of the signal. For the 5° — >■ D^*^^K*^ 
the expected signal to background ratio is worse and we require | cos(6'i^.)| > 0.5. 

We reconstruct D° mesons in the decay channels: K'n'^, K'ti^tt'ti^ and K'n^n^, 
using requirements that the invariant mass be within 20 MeV/c^, 15 MeV/c^ and 25 MeV/c^ 
(~ 3(t) of the nominal mass, respectively. In each channel we further define a mass 
sideband region, with a width twice that of the signal region and location within 0.1 GeV/c^ 
on the nominal mass. For the vr*^ from the K^tt^tt^ decay, we require that its 

momentum in the CM frame be greater than 0.4 GeV/c in order to reduce combinatorial 
background. D*^ mesons are reconstructed in the D*^ —>■ D^n^ decay mode. The mass 
difference between D*^ and candidates is required to be within 4 MeV/c^ of the expected 
value (~ 4(7). 

We combine D^*^^ candidates with Kg or candidates to form B mesons. Signal 
events are identified by their CM energy difference, AE = [Y^^Ei) — E^, and the beam 
constrained mass, Mbc = \/ E^ — i^iPiYi where E\^ is the beam energy and pi and Ei are 
the momenta and energies of the B meson decay products in the CM frame. We select 
events with Mbc > 5.2 GeV/c^ and \AE\ < 0.2 GeV, and define a B signal region of 
5.272 GeV/c^ < Mbc < 5.288 GeV/c^ and \AE\ < 0.03 GeV. In the rare cases where there 
is more than one candidate in an event, the candidate with the L)*^*)" and K^*")^ masses 
closest to their nominal values is chosen. We use Monte Carlo (MC) simulation to model 
the response of the detector and determine the efficiency [lol |. 

To suppress the large combinatorial background dominated by the two-jet-like e+e" — > qq 
continuum process, variables that characterize the event topology are used. We require 
I cos6'thr| < 0.80, where 6'thr is the angle between the thrust axis of the B candidate and that 
of the rest of the event. This requirement eliminates 77% of the continuum background and 
retains 78% of the signal. We also construct a Fisher discriminant, JF, which is based on the 
production angle of the B candidate, the angle of the B candidate thrust axis with respect 
to the beam axis, and nine parameters that characterize the momentum fiow in the event 



relative to the B candidate thrust axis in the CM frame [ll|. We impose a requirement on 
T that rejects 67% of the remaining continuum background and retains 83% of the signal. 
For the suppressed B'^ — > D^*)^K*'^ modes we make a stricter requirement, which has 45% 
efficiency and 93% background rejection. 

Among other B decays, the most serious background comes from B^ D^tt'^, 
D- K^*^^K-n'^, K'-*^^K-n-n+ and B^ , K^*^^n-, i?W°7r-7r°, 

TT TT 7r+. These decays produce the same final state as the 5° D^*^^K^*^^ signal, 
and their product branching fractions are up to ten times higher than those expected for 
the signal. To suppress this type of background, we exclude candidates if the invariant mass 
of the combinations listed above is consistent with the D~ hypothesis within 25 MeV/c^ 
(~ 3(t). The B^ D*^K^, D*^ D^n^ decay can also produce the same final state as 
the B^ — D^K*^ decay; however this decay is kinematically separated from the signal - the 
invariant mass selection criteria for K*^ candidates completely eliminates this background. 
Another potential BB background comes from the D^*^^p^ decay channel with one 

pion from the p° decay misidentified as a kaon. The reconstructed K~tt~^ invariant mass 
spectra for these events overlap with the K*'^ signal mass region, while their AE distribu- 
tion is shifted by about 70 MeV/c^. We study this background using MC simulation. The 
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FIG. 1: AE (left) and Mbc (right) distributions for the 5° ^ i:)OA'(*)o candidates. Points with 
errors represent the experimental data, hatched histograms show the mass sidebands and curves 
are the results of the fits. 



contribution to the 5° — *• signal region is found to be less than 0.2 events. We 

examined the possibility that other B meson decay modes might produce backgrounds that 
peak in the signal region by studying a large MC sample of generic BB events. No other 
peaking backgrounds were found. 

The AE and Mbc distributions for B^ /}0^(*)o candidates are presented in Fig. ^ 
where all three decay modes are combined. Each distribution is made for events from 
the signal region of the other parameter. The hatched histograms in Fig. ^ are the dis- 
tributions for events in the mass sideband scaled according to the D° mass selection. 
The sideband shape replicates the background shape well, confirming that the background 
is mainly combinatorial in nature. Clear signals are observed for the D^K^ and D^K*^ final 
states. As a cross check, we also study the candidates' invariant mass and flight distance 
distributions and K*^ candidates' invariant mass and helicity distributions for these decays. 
The distributions mentioned above are shown in Fig. |2l where the points with error bars are 
the results of fits to the AE spectra for data in the corresponding bin, and the histograms 
are distributions from signal MC. All distributions are consistent with the MC expectations. 

For each decay mode, the AE distribution is fitted with a Gaussian for signal and 
a linear function for background. The Gaussian mean value and width are fixed to the 
values from MC simulations of the signal. The region AE < —0.1 GeV is excluded from 
the fit to avoid contributions from other B decays, such as i? — > D^*''^K^*''^{ti) where (vr) 
denotes a possible additional pion. For the Mbc distribution fit we use the sum of a signal 
Gaussian and an empirical background function with a kinematic threshold with a 
shape parameter fixed from the analysis of the off-resonance data. For the calculation of 
branching fractions, we use the signal yields determined from fits to the AE distribution. 
This minimizes the possible bias from other B meson decays, which tend to peak in Mbc 
but not in AE. The fit results are presented in Table U where the listed efficiencies include 
intermediate branching fractions. The statistical significance of the signal quoted in Table H] 
is defined as a/— 2 ln(£o / J^max) , where £max and Cq denote the maximum likehhood with 
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FIG. 2: From left to right: Distributions of invariant mass and flight distance for the Kg candidates' 
in — > D^K^ channel, invariant mass and helicity distributions for the K*'^ candidates' in 
^0 ^ £)0^*o channel. 
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FIG. 3: /^E distributions for the 5° ^ D*^K^*)^ and 5° ^ D^*)^K*^ candidates. Points with 
error bars represent the experimental data, hatched histograms show the mass sidebands and 
curves show the results of the fits. 

the nominal signal yield and the signal yield fixed at zero, respectively. 

For the final result we use a simultaneous fit to the /S.E distributions for the three decay 
channels taking into account the corresponding detection efficiencies. The normalization of 
the background in each sub-mode is allowed to float while the signal yields are required 
to satisfy the constraint Ni = N^g ■ B{B^ — > • , where the branching fraction 

B{B^ is a fit parameter; A''^^ is the number of BB pairs and Si is the efficiency, 

which includes all intermediate branching fractions. 

The statistical significance for S° D*^K^ is 3.2cr. This is the first evidence for this 



TABLE I: Fit results, efficiencies, branching fractions and statistical significances for 
^(*)o^(*)o decays. 



Mode 




AE yield 


EfBciency (IQ-^) 


B (10-5) 


Significance 






78.1 ±14.1 


8.02 


3.72 ±0.65 ±0.37 


6.6a 




D°K*° 


77.7 ±14.6 


9.98 


3.08 ±0.56 ±0.31 


6.3a 




D*°K° 




1.97 


3.18li:i^±0.32 


3.2a 








2.54 


2.Mtlit{< 4-8) 90% CL 


2.1a 




£,0^*0 




6.52 


< 0.51 90% CL 








0+2.7 
o.o_2.i 


1.72 


< 1.9 90% CL 





7 



decay. We do not observe significant signals for tlie —>■ D*^K*^ and tlie /}(*)07^*o 
decay modes and set 90% CL upper limits for them. Figure 01 shows the distributions for 
S° D*^K^*'^^ and 5° Di*)^K*^ candidates. The upper limit N is calculated from the 
relation C{n)dn = 0.9 Jq°° C{n)dn, where C{n) is the maximum likelihood with the signal 
yield equal to n. We take into account the systematic uncertainties in these calculations by 
reducing the detection efficiency by one standard deviation. Using the measured branching 
fraction for B^ D^K*^ decay and upper limit for B^ D^K*^ ^ we set an upper limit on 
the amplitude ratio: Rdok-o < 0.39. 

As a check, we apply a similar procedure to the decay chains with the same final states: 
5° D+[K^K+]n^ and B^ D*+[D^7i+]K-. The estimated branching fractions are 
consistent with the world average values jl3i] . 

The following sources of systematic errors are found to be significant: tracking efficiency 
(1-2% per track), kaon identification efficiency (1%), vr^ efficiency (6%), Kg reconstruction 
efficiency (6%), efficiency for slow pions from D^-K^ decays (8%), D^*'^^ branching 

fraction uncertainties (2% - 6%), signal and background shape parameterization (4%) and 
MC statistics (2% - 3%). The tracking efficiency error is estimated using rj decays to 77 
and 7r"'"7r~7r°. The kaon identification uncertainty is determined from D*^ — > D^vr^, 
K^TT^ decays. The vr*^ reconstruction uncertainty is obtained using decays to K^n^ 
and K^-K^-K^ . We assume equal production rates for B^B^ and B^B^ pairs and do not 
include the uncertainty related to this assumption in the total systematic error. The overall 
systematic uncertainty is found to be 10% for 5° ^ D°K^*^° and 13% for 5° ^ D*°K^*^^. 

In summary, we report improved measurements of 5° branching fractions 

B{B° D^K^) = (3.72 ± 0.65 ± 0.37) x 10"^ and B{B° D°K*°) = (3.08 ± 0.56 ± 
0.31) X 10~^. Note that we ignore the possible contribution of B^ —>■ D^K^ to the former 
result, since we do not distinguish between and K^. We also report the first evidence, 
with 3.2(7 significance for B^ — >■ D*^K^, with a branching fraction of B{B^ D*^K^) = 
{S.lStlil ± 0.32) X 10"^ No significant signal is observed in the B° jj*ok*o final state. 
The corresponding upper limit at the 90% CL is i3(5° d*^K*^) < 2Mt\ll x 10^^- We 
also set the 90% CL upper limits for the Vub suppressed 5° — > Z)(*)o^*° decays: B{B'^ 
D^K*^) < 0.51 X 10-5 and BiB^ ^*o^*0) ^ ^ g ^ ^q-5 r^j^g results obtained supersede 
our previous measurements ; these results are consistent with, and more precise than, the 
previous measurements. 

We thank the KEKB group for the excellent operation of the accelerator, the KEK 
Cryogenics group for the efficient operation of the solenoid, and the KEK computer group 
and the National Institute of Informatics for valuable computing and Super-SINET network 
support. We acknowledge support from the Ministry of Education, Culture, Sports, Science, 
and Technology of Japan and the Japan Society for the Promotion of Science; the Australian 
Research Council and the Australian Department of Education, Science and Training; the 
National Science Foundation of China under contract No. 10175071; the Department of 
Science and Technology of India; the BK21 program of the Ministry of Education of Korea 
and the CHEP SRC program of the Korea Science and Engineering Foundation; the Polish 
State Committee for Scientific Research under contract No. 2P03B 01324; the Ministry of 
Science and Technology of the Russian Federation; the Ministry of Education, Science and 
Sport of the Republic of Slovenia; the Swiss National Science Foundation; the National 
Science Council and the Ministry of Education of Taiwan; and the U.S. Department of 
Energy. 



8 



* on leave from Nova Goriea Polytechnic, Nova Gorica 
[1] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973). 
[2] I. Dunietz, Phys. Lett. B 270, 75 (1991). 

J.H. Jang, P.Ko, Phys. Rev. D 58, 111302 (1998). 
[3] M. Gronau, D. London, Phys. Lett. B 253, 483 (1991). 

B. Kayser, D. London, Phys. Rev. D 61, 116013 (2000). 
[4] Belle Collaboration, P. Krokovny et al, Phys. Rev. Lett. 90, 141802 (2003). 
[5] M.Gronau Phys. Lett. B 557, 198, (2003), 

R. Fleischer, Phys. Lett. B 562, 234 (2003). 
[6] The inclusion of charge conjugate states is implicit throughout this report. Note that we do 

not distinguish between final states with and K^. 
[7] Belle Collaboration, A. Abashian et al, Nucl. Inst, and Meth. A 479, 117 (2002). 
[8] S. Kurokawa and E. Kikutani, Nucl. Instr. and Meth. A 499, 1 (2003). 
[9] Y. Ushiroda (Belle SVD2 Group), Nucl. Instr. and Meth. A 511 6 (2003). 
[10] R. Brun et al, GEANT 3.21, CERN DD/EE/84-1, 1984. 
[11] CLEO Collaboration, D.M. Asner et al, Phys. Rev. D 53, 1039 (1996). 
[12] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B 241, 278 (1990). 
[13] S. Eidelman et al (Particle Data Group), Phys. Lett. B 592, 1 (2004). 



9 



